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Abstract:  This  paper  adresses  the  problem  of  determining  upper  and  lower  bounds  for 
the  effectivity  index  on  the  a-posteriori  estimate  of  the  error  in  the  finite  element  method. 
These  bounds  axe  given  explicitly  for  a  certain  concrete  estimator  for  linear  elements  and 
unstructured  triangular  meshes.  They  depend  strongly  on  the  geometry  of  the  triangles 
and  (relatively  weakly)  on  the  smoothness  of  the  solution.  An  example  shows  that  the 
bounds  are  not  over  pessimistic.  In  [4]  detailed  numerical  experimentation  is  given. 
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1.  Introduction.  Since  the  first  papers  by  Babuska  and  Rheinboldt  on  the  a-poste- 
riori  estimation  of  the  errors  in  the  finite  element  method  [5,6],  this  subject  became  an 
increasingly  important  aspect  of  the  application  of  this  method.  During  the  last  years  sev¬ 
eral  codes  including  different  estimators  have  been  developped  [14,23,25,26,28]  and  nowa¬ 
days  there  are  many  different  estimators  in  use  for  a  given  problem  (see,  for  instance, 
[12,13,21,24]  and  references  there  in). 

A  standard  measure  of  the  quality  of  an  estimator  is  the  so  called  effectivity  index 


eff  = 


estimated  error 
true  error 


For  a  given  problem  an  estimator  is  said  to  be  equivalent  to  the  error  if  the  effectivity 
index  is  bounded  below  and  above  by  two  strictly  positive  constants  independently  of  the 
meshsize: 


c<eff<C; 


these  constants  may  depend  on  the  class  of  functions  under  consideration.  (Here  and 
thereafter,  c  and  C  will  denote  constants  not  necessarily  the  same  at  each  occurrence,  but 
always  independent  of  the  meshsize). 

A  property  that  has  been  considered  highly  relevant  to  measure  the  potential  quality 
of  an  estimator  is  the  so  called  asymptotic  exactness.  Roughly  speaking,  an  estimator  is 
asymptotically  exact  for  a  particular  problem  if  its  effectivity  index  converges  to  one  when 
the  meshsize  aproaches  to  zero. 

In  the  one  dimensional  case  Babuska  and  Rheinboldt  [7,8]  made  a  complete  analysis 
of  asymptotically  exact  error  estimators.  For  two  dimensional  elliptic  problems,  several 
estimators  have  been  proved  to  be  asymptotically  exact  when  used  on  almost  uniform 
patches  of  rectangular  or  triangular  meshes,  provided  the  solution  of  the  problem  is  smooth 
enough  [2,11,17,18,19]. 

In  particular,  for  linear  triangular  elements,  some  well  known  local  estimators  like 
Bank-Weiser’s  [15]  and  Zienkiewicz-Zhu’s  [31]  are  asymptotically  exact  on  uniform  meshes 
as  that  in  Figure  1.1. a  but  not  on  other  rather  uniform  meshes  as  those  in  Figures  1.1. b 
and  1.1. c.  (See  [19]  for  Bank-Weiser’s  estimator  and  [9]  for  Zienkiewicz-Zhu’s). 
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(a)  (b)  (c) 

Figure  1.1 

A-posteriori  error  indicators  (i.e.:  estimators  per  element)  are  employed  in  adaptive 
processes  to  identify  those  portions  of  the  mesh  with  bigger  errors  in  order  to  generate  a 
new  refined  mesh.  Usually,  the  meshes  generated  by  these  adaptive  processes  axe  regular 
(in  the  sense  of  a  minimal  angle  condition)  but  not  uniform  as  in  Figure  1.1. a.  Very  likely, 
Jill  the  used  estimators  are  not  asymptotically  exact  on  the  meshes  that  axe  adaptively  con¬ 
structed.  However,  the  estimators  actually  in  use  are  equivalent  to  the  error  for  any  regular 
family  of  meshes  with  bounds  on  the  effectivity  index  depending  only  on  the  regularity  of 
the  mesh.  Anyway,  in  no  case  these  bounds  axe  known  explicitly.  To  increase  the  accuracy 
of  the  indicators  and  estimators,  various  correction  factors  derived  by  computational  tests 
are  used. 

In  this  paper  we  shall  analyze  a  particular  estimator  based  on  Babuska-Miller’s  [3]; 
(this  type  of  estimator  is  used,  for  instance,  in  [25]).  We  shall  prove  again  the  equivalence 
of  this  estimator  for  the  Laplace  equation,  but  in  such  a  way  that  it  will  be  able  to  compute 
asymptotic  bounds  of  its  effectivity  index  in  terms  of  the  geometry  of  the  mesh  and  on 
the  smoothness  of  the  solution.  We  shall  show  that  these  bounds  axe  sharp  and  that  their 
dependence  on  the  geometry  of  the  mesh  is  optimal.  Finally  we  shall  present  similar  results 
for  the  elasticity  problem. 

2.  The  error  estimator.  Let  us  consider  as  our  first  model  problem  the  Laplace 
equation  with  mixed  boundary  conditions.  Let  12  be  a  bounded  polygon  in  R2  and  let 
its  boundary  df2  be  split  into  two  parts  I’d  and  Tn  (I'd  of  positive  length).  Let  u  be  the 
solution  of  the  problem 


(2.1) 


(  —Au  =  f  ,  in  0,  , 
<  u  =  0  ,  on  Td  , 


t  =3 


on  rn , 


4 


where  n  is  the  unit  outer  normal  vector  to  dQ,  f  €  L2(Q)  and  g  €  L2(rn). 

We  shall  use  the  standard  notation  for  Sobolev  spaces  Hm(D),  their  norms  |j  •  ||m  D 
and  seminorms  |  •  |m  D.  Let  -Hpd(f2)  :=  {u  €  H1^)  :  u|pd  =  0}.  |  •  ]j  Q  is  a  norm  on  that 
space;  it  is  the  energy  norm  of  this  model  problem. 

Let  {Tft}  be  a  regular  family  of  triangulations  of  Q  (i.e.:  the  minimal  angle  of  all  the 
triangles  is  bounded  below  by  a  positive  constant,  the  same  for  all  the  meshes);  as  usual  h 
stands  for  the  maximal  meshsize  and  we  assume  that,  when  the  edge  of  a  triangle  intersects 
5f2,  it  is  completely  contained  either  in  I'd  or  in  Fn.  The  meshes  axe  not  assumed  to  be 
quasiuniform. 

Let  u/,  €  14  :=  {v  €  H^d{?l)  :  v\T  €  V\ (T),  VT  €  Th)  be  the  piecewise  linear  finite 
element  approximate  solution  of  problem  (2.1).  ( Vm(T )  denotes  the  set  of  polynomial 
functions  defined  on  T  of  degree  not  greater  than  m).  Let  e  :=  u  —  u/,  denote  the  error  of 
this  approximation. 

Integrating  by  parts  we  obtain  for  any  v  €  H^d(Q)  : 

f  Ve  •  Vu  =  f  Vu  •  Vu  -  V  f  Vufc  •  Vu  =  f  fv+  f  gv  -  V  /  , 

Jn  Jn  Jt  J n  J rn  Jar  dnT 


where  for  each  triangle  T,  nT  is  its  unit  outer  normal  vector. 

Let  us  call  T;  the  union  of  all  the  interior  edges  of  the  triangulation  Th.  For  each  edge 
£  C  T;  let  us  choose  an  arbitrary  normal  direction  n  and  denote  the  two  triangles  sharing 
this  edge  T\n  and  Tou t,  where  the  normal  n  is  outwards  T\n.  Let 


duh 

dn 


:=  v(u'*lrou,)-n-  V(ufc|Tin)-n 


denote  the  jump  of  accross  the  edge  £;  this  value  is  independent  of  the  choice  of  n. 
With  this  notation  we  may  now  write  the  so  called  residual  equation: 

(2-2)  Lv'  v°-Jjv+LXg-d^)v+£jMb' Wv€HKin)- 


This  equation  relates  the  error  e  with  the  interior  residuals  f\T  =  —  [A(u|r)  —  A(uft|r)]  , 
the  boundary  residual  (g  —  ^*-)  and  the  jumps  of  the  gradient  of  the  finite  element  ap¬ 
proximation  [^-]r  Several  estimators  have  been  obtained  by  approximating  the  error  as 
the  solution  of  this  equation  [3,15,16,29].  The  estimator  that  we  shall  consider  is  a  slight 
variation  of  Babuska-Miller’s  [3]  that  Verfurth  describes  for  the  Stokes  problem  [29]. 

For  any  triangle  T  €  Th,  let  Er  be  the  set  of  its  three  edges  and  let 
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be  the  £2(T)-projection  of  /  onto  the  constants.  For  any  edge  t  <=  rn,  let 


n,g:  nl! 


be  the  X2(^)-projection  of  g  onto  the  constants.  For  each  edge  £  of  the  triangulation,  let 

lo,  if/crd. 


(2.3) 


J,  ■■=  { 


We  define  as  an  estimator  of  the  local  energy  error  |e[j  T, 


(2.4) 


TJt  ■ — 


m2(nr/)!  +  i  e 

(€Et 


Although  we  deal  with  the  Laplace  model  problem,  this  approach  is  valid  for  any 
divergence  type  operator  with  piecewise  constant  coefficients  if  the  meshes  axe  such  that 
the  interfaces  of  the  coefficients  coincide  with  boundaries  of  elements. 

3.  Equivalence  between  the  error  and  the  estimator.  The  ideas  of  Verfiirth 
[30]  can  be  directly  applied  to  our  simpler  model  problem  to  prove  the  following  theorem 
without  any  further  assumption  on  the  mesh  and  for  any  problem  (2.1)  with  solution 
u  €  H'iSl). 

THEOREM  3.1.  With  the  definitions  and  assumptions  of  Section  2,  there  exist  two 
positive  constants  C  and  C'  only  depending  on  the  regularity  of  the  mesh  such  that 


(3.1) 

leli,n  —  C 

E 

U?  +  m  ll/-nr/H*r+  E  Mlls-n.jIlL] 

T€T* 

\  *COTnr„)  / 

and 

(3-2) 


TJr<C' 


K,f+  E  inil/-nr,/||02r,  + 


\T’CT 


r 


E  mi  Is-n.jllb 

,/c(arnrn) 


where  T  :=  IJ{T'  €  7*  :  T  and  T'  have  a  common  edge}  . 

Proof.  The  proof  will  not  be  given  here  because  it  is  essentially  identical  to  that  in 
[30].  0 
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These  bounds  show  that  whenever  the  data  /  and  g  axe  locally  smooth,  if  the  error  is 
properly  0(hs)  (0  <  s  <  1),  then  the  estimator 

(3.3)  In  ==  (  £  ’ll) 

\Terk  / 

is  globally  equivalent  to  the  error.  In  fact,  we  have  the  following  theorem. 

THEOREM  3.2.  In  addition  to  the  assumptions  of  Section  2,  let  us  assume  that  there 
exists  a  triangulation  T  such  that 

(3.4)  f\T£Hl(T),  VT  €  T 
andVTeT  :  dT  n  T„  #  0 

(3.5)  g\t  6  Hl{t)  ,  W€3Tn  Tn  ; 

let  us  also  assume  that  all  the  triangulations  Th  are  refinements  of  T. 

If  there  exist  constants  C*  >  0  and  s  €  (0, 1]  not  depending  on  h  such  that 

(3.6)  |e|j  >  C*ha , 
then  there  exist  two  positive  constants  c  and  C  such  that 

(3-7)  crjn  <  |e|j >n  <  Cr)n  . 

Proof.  By  using  (3.1-5),  the  regularity  of  the  meshes  and  the  standard  aproximation 
properties  of  the  projections  IIT/  and  II,  <7,  we  may  write 


since  each  T  is  the  union  of  at  most  4  triangles  of  7/,.  Hence,  by  using  (3.6)  the  theorem 
is  proved.  □ 
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Remark  3.1.  The  assumptions  about  the  existence  of  T  is  made  only  to  cover  those 
cases  where  /  and  g  are  piecewise  smooth  and  the  meshes  sire  such  that  the  interfaces  of  the 
data  coincide  with  boundaries  of  the  elements.  On  the  other  hand,  these  local  smoothness 
assumptions  (3.4)  and  (3.5)  can  be  weakened;  in  fact,  if  f\T  €  Hf(T)  and  g\(  €  H $  +<(£) 
for  some  e  >  0,  then  the  conclusion  of  Theorem  3.2  and  all  what  follows  are  valid.  □ 


Remark  3.2.  The  error  is  always  properly  0(ha )  (i.e.:  assumption  (3.6)  is  valid) 
except  for  trivial  cases  (see  [3]).  □ 

We  shall  now  describe  a  variation  of  Verfurth’s  proof  of  Theorem  3.1  that  will  give 
computable  asymptotic  approximations  of  the  constants  c  and  C  in  the  equivalence  (3.7) 
(assuming  slightly  stringent  hypothesis  for  the  upper  bound).  In  the  following  sections,  we 
shall  show  that  the  constants  obtained  in  this  way  are  almost  achievable. 

Let  Vp  :=  {v  €  -Hpd(f2)  :  v\T  €  VP(T),  VT  6  T/,}  (in  particular  Vi  =  V/,)  and,  forp  >  2, 
let  up  €  Vp  be  the  finite  element  approximate  solution  of  problem  (2.1)  in  this  space.  Let 
tp  :=  Up  —  Uh ,  then 

(3.8)  ieli,n  ^  lu  —  upli,n  +  lepli,n  • 


It  is  known  [10]  that,  if  u  €  Hl+((Q.)  for  some  e  >  0  and  if  the  family  of  meshes  is 
quasiuniform,  then 


(3.9) 


i«-«Pilin  <  C7imin<p’ey 


u 


l  +  e,0 


•> 


with  a  constant  C  independent  of  u,  h  and  p.  Therefore,  if  the  solution  is  smooth  enough, 
say  u  6  H1+t(Q)  for  some  e  >  1,  |u  —  liplj  n  is  asymptotically  negligible  with  respect  to 
the  error  lelj  n  for  any  p  >  2.  Instead,  if  the  solution  u  €  for  some  e  €  (0, 1],  the 

error  |e|j  n  is  expected  to  be  0(hf)  and,  in  this  case,  |u  —  up|j  n  will  be  negligible  with 
respect  to  le^  n  only  for  p  big  enough.  In  any  case,  even  for  h  small  or  for  p  big  enough 
(or  both  together),  the  term  |it  —  up\1  n  can  be  neglected  in  (3.8).  Therefore,  it  is  enough 
to  bound  |epjj  Now, 

(3.10)  |ep|j  n  =  /  Vep  •  Vep  =  /  V(up-u)-Vep+  /  Ve  •  Vep  =  /  Ve  •  Vep  , 

’  Jn  J  n  Jo,  Jo 

where  we  have  used  that  ep  €  Vp. 

For  any  continuous  function  v  defined  on  SI,  let  v1  denote  its  Lagrange  piecewise  linear 
interpolant  on  the  mesh  T*.  Since  €  V*,  then  from  (3.10)  we  have 


and  by  using  the  residual  equation  (2.2)  we  may  write 


where 

«(»):=  £  f  (f-  nr/)t>  +  £  f(g-n,g)v. 

T€Th  Jt  ec rn  Jt 

Using  the  definition  (2.4)  of  the  local  estimator  t]t  and  Cauchy-Schwarz  inequality  we 
have: 


(3.11)  |ep|j  n  <  Vt 

T€Th 


\T\  St'6*  e'r)  +  2  £  |  \t\  J,{e '  4 


+  KeP  -  ej)| 


This  last  expression  allows  us  to  prove  the  following  theorem. 
THEOREM  3.3.  Under  the  assumptions  of  Theorem  3.2, 


eli,n  - 


£  nfril 


(3.12) 


TZTh 

+  C 


+  |u-up|in 


(  5Z  I/I It)  +fcf 

rer  /  \iec 


in 

2 


where 


(3.13)  (C')2:=  sup 

v€Vp\V0 


ITT  JAv~vI)  +12Zt<=ET 


ii,r 


Proof.  According  to  (3.11)  and  the  definition  of  C£: 


lepli,n  ^  E  +  I s(ep  ~  ei)|  £ 

T€Th 


E  (CPr) 

L  rerfc 


2  2 


lepli,n  +  l^(ep  ep)|  • 
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Proceeding  as  in  Theorem  3.2  we  prove  that 


i««,  -  «;)i  <  c  +*|(Ewh') 

[  Vrer  /  \tec  / 

So,  by  using  (3.8)  we  conclude  the  theorem.  Q 

The  constants  C£  in  this  theorem  depend  on  the  degree  p  used  to  make  \u  —  up\l  u 
negligible  in  (3.12).  However  the  next  theorem  shows  that  this  dependence  is  very  weak. 

THEOREM  3.4.  Let  CP  be  defined  by  (3.13);  there  exists  a  constant  CT  only  depending 
on  the  shape  of  the  triangle  T  such  that  Vp  >  2 

CP  <  CT  log^p 

Proof.  Let  T  {(ar,y)  :  x  >  0,  y  >  0,  and  x  +  y  <  1}.  For  any  polynomial  v  of 
degree  p  >  2 

H*?Ht«(f)  -  Clo&P  Flli,f 

with  C  independent  of  v  and  p;  (this  is  an  inmediate  consequence  of  Theorem  6.2  in  [1]). 
Since  ( v  —v1)  vanishes  for  v  G  F0,  then 

11^  —  ^/H£,oo(f?)  -  CloS ’ 

and  so,  for  any  v  G  VP(T),  by  changing  coordinates  to  the  triangle  T  we  obtain 

lb-  v/IL~(T)  ^  ct  losv 

with  a  constant  CT  only  depending  on  the  shape  of  the  triangle.  Using  this  inequality  in 
the  definition  (3.13)  of  C£  we  conclude  the  theorem.  0 

In  the  following  section  we  shall  compute  the  costants  C£  for  different  values  of  p  and 
we  shall  analyze  their  dependence  on  the  shape  of  the  triangle  T.  On  the  other  hand,  for 
the  lower  bound  in  (3.7)  we  have  the  following  theorem. 

THEOREM  3.5.  For  each  mesh  Th ,  let  w  G  H(d(fl)  (eventually  depending  on  the  mesh) 
be  such  that  for  all  the  triangles  T  €  Th, 

(3.14)  /  (nT/)u>  =  |T|2(nT/)2 , 

JT 


(3.15) 


E  =  E 

c  r-  C  J  '  tr- 


2 

t 
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and 


(3.16) 


3c;  >0  :  H,,T<c;,r , 


where  C'T  may  depend  on  the  shape  of  the  triangle  but  not  on  its  size  hT . 
Then,  under  the  assumptions  of  Theorem  3.2, 

(3.17)  <  (sup  CTJ  •  |e|,,n  +  C  E  +  hi  (E  l?lb) 


Proof.  By  using  (3.14)  and  (3.15)  in  the  definition  of  r)T ,  the  residual  equation  (2.2) 
and  the  definition  of  6,  we  have 


Etf-E 

T€Tk  T£Tk 


and  hence, 


f  (IlTf)w  +  53  I  Jew  =  [  Ve  •  Vw  - 

2  (tEr Je 


53  +  l*(«OI 

Ter* 


6(w ) 


Now 


E  /(/-  nr/)w  =  Y,  /(/-nr/)(»-nT») 

T€T*  '/t  T€T*  Jt 


<c£|T| 

T£Tk 


and 


53  /(^~n^)w  =  53 


<  c  53  I*I*Mi,<Mi,7>  > 

/crn 


where  T(  is  the  triangle  in  7),  such  that  i  C  dTt .  Therefore,  by  using  (3.16),  we  have 

$ 


53  -  {  leli,n  +  C 

T€T» 


<  (  sup  C '  }  < 

\T€Tfc  / 


lcli,n  +  c 


h2  £l/!?,r  +^3  EWm 

\T€T  /  \*€£ 

*2fEi/Hr)  +^(ei<< 

\T€T  /  \*€C 


>  W 


i.n 


r 


E*i 

<T€Th 


and  hence  we  obtain  (3.17). 


□ 


In  the  next  section  we  shall  exhibit  functions  w  satisfying  the  hypothesis  of  Theorem 


3.5  and  we  shall  show  how  to  calculate  the  constant  C'  . 

T 
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4.  Computation  of  the  bounds.  In  order  to  compute  the  constants  C£  of  Theorem 
3.3,  let  zq  6  'Pp(T)  be  the  solution  of  the  weak  finite  dimensional  problem 

(4.1)  JtVz0.Vv=JL  Vv<iVp(T) 

and,  for  i=l ,2,3,  let  z,  €  VP(T)  be  the  solution  of 


<4'2)  JTVzi-Vv  =  7k\l{'’-v')’  VveVr(T)’ 

where  £,,  i  =  1,2,3,  are  the  three  edges  of  T.  We  may  write 


(CP)2  =  sup 
v€-Pp\Vo 


Elo 


ll,T 


=  sup 

v€2\P o 


Elo  Utv*‘-V”Y 

\v\l,T 


where  Z  is  the  subspace  of  VP(T )  spanned  by  {2i}^_0. 

3  2 

For  v  =  2Zj=o  u«2‘  €  Z  we  may  write  lulj  T  =  vlCv,  where  v  :=  (v0,...,W3)  and 
C  €  R4x4  is  the  symmetric  matrix  of  entries  Cij  :=  fT  Vzj  •  Vzj,  i,j  =  0, ...  ,3.  On  the 
other  hand  ]T^_0  (Jt^z‘  ‘  ^ v Y  =  vtC2v.  Therefore, 


(c  py 


vtC2-' 


is  the  spectral  ratio  of  C. 

So,  to  compute  the  constants  for  any  degree  p  >  2  and  any  triangle  T,  we  only  need  the 
solutions  z,  of  problems  (4.1)  and  (4.2).  These  functions  are  the  p-degree  finite  element 
solutions  of  elementary  elliptic  problems  on  the  triangle  T  with  a  mesh  consisting  of  this 
only  triangle;  they  have  been  computed  by  using  the  code  PROBE  [27].  Our  computations 
show  that  for  any  triangle  T  and  for  any  degree  p  =  2, 3, . . . ,  8, 

(4.3)  0.548  log^p  sin'^^j  <  C£  <  0.813  log^p  sin~^^y-j  , 

where  aT  is  the  minimum  angle  of  T.  These  constants  CPT  also  depend  on  the  other  angles 
of  T;  however  this  dependence  is  very  weak.  In  fact,  the  estimate  (4.3)  is  valid  for  all  the 
triangles  with  minimum  angle  aT,  independently  of  the  size  of  the  other  angles. 

From  Theorem  3.4  we  know  that  for  any  fixed  triangle  the  constants  Cp  jure  bounded 
above  by  log ^ p;  our  computations  show  that,  actually,  they  are  almost  proportioned  to 
log^p.  On  the  other  hand,  for  a  fixed  degree  p  >  2,  the  constants  depend  on  the  geometry; 
they  essentially  depend  on  the  minimum  angle  and  in  fact  they  deteriorate  when  this  angle 
is  very  small,  but  the  square  roots  in  (4.3)  makes  this  dependence  to  be  weak. 
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Now,  we  shall  describe  how  to  compute  the  constants  C^.  of  Theorem  3.5.  To  this 
goal  we  need  a  function  w  €  satisfying  (3.14),  (3.15)  and  (3.16)  with  constants 

C'T  as  small  as  possible.  We  define  this  function  w  in  each  triangle  but  in  such  a  way  that 
it  satisfies  the  required  global  smoothness.  For  any  edge  i  of  the  triangulation  we  choose 
a  continuous  function  tj>t  vanishing  at  both  ends  of  the  edge  and  such  that  its  average 
ct  :=  jlj  ^  0.  To  guarantee  that  w  €  we  consider  only  those  functions  w 

whose  restrictions  to  t  are  a  multiple  of  ij)l  satisfying  (3.15);  therefore  tu^  =  • 

We  shall  introduce  some  notation  in  order  to  define  w  in  the  interior  of  each  triangle 
T.  Let  r  :=  (IIT/,  Jly,  J*2,  Je3)  G  R4  and  D  :=  diag(|T|,  -^=,  -^);  then  r/2  =  rlD2r. 

Let 


WJ 


:=j  weH\T)  :  Jw  =  \T\2(UTf)  and  w\ti  =  ,  *  =  1, 2,  sj  ; 


Wj  is  an  affine  subspace  of  Hl(T)  parallel  to  the  subspace 

Kj  :=  |iu  £  HX{T)  :  J  w  =  0  and  w|<(.  =0,  i  =  1,2, 3^ 
Let  wJ  €  Wj  be  such  that 


(4.4) 


min 

u 


then  wj  satisfies  (3.14),  (3.15)  and  (3.16)  with  a  constant 


sup 

\r#0 


rlD2r  J 


1 

2 


To  compute  this  constant  we  need  to  calculate  /j’|Vu;^|2  for  any  r  €  R4.  Let  us 
remark  that  (4.4)  holds  if  and  only  if  €  Wj  satisfies  fT  WwJ  •  Vu;  =  0,  Vw  €  Kj. 

Let  wo  be  the  solution  of  the  Dirichlet  problem 


(4.6) 


j  —  Au>o  =  1  ,  in  T  , 

1  wqIqt  =  0  , 


then 


Vu> o  •  Vio 


dw0 

~dn 


For  i  =  1,2,3,  let  w,  be  the  solution  of 


Vw  €  Kj . 


(4.7) 


J  —Awt  =  0  ,  in  T  , 

1  Wi\l{  —  li>ti  ,  ^ildTV;  =  0  1 
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(we  have  used  that,  because  of  (4.6),  fT  |Vtuo|2  =  fTwo). 
Finally,  because  of  (4.7),  fTVwi  •  Vw0  =  0,  and  so 


which  is  a  quadratic  form  on  r.  Therefore,  the  computation  of  the  constant  C'T  by  means  of 
(4.5)  reduces  to  a  simple  eigenvalue  problem  which  can  be  easily  solved  once  the  solutions 
i vi  of  the  Dirichlet  problems  (4.6)  and  (4.7)  are  known.  In  our  computations  we  have  also 
used  the  code  PROBE  to  solve  numerically  these  problems. 

The  function  w  €  obtained  by  patching  together  all  the  wj  for  T  €  Tfc,  gives 

the  best  possible  constants  for  each  triangle  for  a  given  choice  of  the  edge  functions  rpt. 
After  some  experimentation  we  choose  ipt  as  quadratic  functions  vanishing  at  both  ends 
of  the  edge.  This  choice  gives  constants  satisfying  for  any  triangle  T : 

(4.8)  3.45 sin"*  (^")  ^  C'T  -  5.85  sin-"*  , 

where  aT  is  the  minimum  angle  of  T.  Once  again,  C^,  is  almost  proportional  to  sin-  2  (~2~) 
and  practically  independent  of  the  size  of  the  othe  angles  of  the  triangle. 

Finally,  by  using  (4.3)  and  (4.8)  and  Theorems  3.3  and  3.5,  we  obtain 


0.171  sin*(|)  rjn  -C  fc2^^l/l?,rj 

(4.9) 

<  0.813  log*p  sin" 2  rjn  +  |u  -  u„|1>n  +  C  h2  ^  ^  |/|?>Tj  +  hi  Iff |J., 
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where  a  is  the  minimum  angle  of  the  mesh  Th ■  The  bounds  (4.9)  cam  be  made  more 
accurate  for  specific  values  of  the  minimal  angle  at  and  of  the  degree  p  >  2;  Table  4.1 
shows  values  of  the  constants  C'a  and  C£  for  the  estimate 


(  \  * 

/  \  * 

C'avn-C 

WEi/iM 

+  '‘l(EWb) 

\T€T  / 

\tec  / 

<  r?n  +  \u  -  +  C 


1 


for  different  values  of  a  and  p. 


a 

C'a 

P  =  2 

p  =  4 

p  —  6 

P  =  8 

7.5° 

2.390 

3.306 

3.660 

3.988 

15.0° 

1.682 

2.341 

2.609 

2.839 

22.5° 

1.363 

1.918 

2.156 

2.343 

30.0° 

0.099 

1.169 

1.670 

1.895 

2.058 

37.5° 

0.108 

1.035 

1.508 

1.727 

1.876 

45.0° 

0.115 

0.939 

1.400 

1.615 

1.757 

52.5° 

0.119 

0.875 

1.334 

1.547 

1.684 

60.0° 

0.121 

0.850 

1.309 

1.522 

1.657 

Table  4.1.  Constants  of  equivalence. 


5.  Sharpness  of  the  bounds.  We  shall  analyize  the  sharpness  of  the  estimates 
obtained  in  the  previous  section  by  considering  a  simple  example.  In  particular,  we  shall 
show  that  the  dependence  of  these  bounds  on  the  geometry  of  the  mesh  is  optimal. 

Let  us  consider  a  particular  case  of  problem  (2.1)  where  Cl  is  a  rectangle  as  in  Figure 
5.1,  rd  consist  of  the  two  vertical  edges  of  Cl  and  Tn  of  the  horizontal  ones;  let  /  be  a 
constant  and  <7  =  0.  The  solution  is  a  quadratic  polynomial  in  x  (and  it  does  not  depend 
on  y).  Let  Th  be  a  family  of  uniform  meshes  like  that  in  Figure  5.1. 
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Since  the  solution  is  quadratic  and  the  Neumann  boundary  conditions  are  zero,  for  any 
of  these  meshes  the  finite  element  approximation  is  exact  at  the  nodes.  Therefore,  it  is 
possible  to  compute  explicitly  the  true  error  and  the  estimator.  The  error  is  the  same  for 
all  the  triangles;  it  only  depends  on  the  meshsize  h  and  on  the  angle  0  which  measures  the 
regularity  of  the  mesh  (see  Fig.  5.1).  For  all  the  elements  disjoint  with  I'd  the  estimator 
is  also  the  same;  for  those  elements  with  an  edge  l  on  the  boundary  Ta,  the  estimator  will 
be  smaller  since,  according  to  (2.3),  the  corresponding  “jump”  Jt  =  0.  However,  since  the 
proportion  of  the  elements  with  an  edge  on  Ta  goes  to  zero  when  the  mesh  is  refined,  the 
global  effectivity  index  is  in  this  case,  asymptotically  equal  to  the  local  one  effr  :=  t  ■ 

An  explicit  computation  gives  eff^  =  18  cot/3.  Let  a  denote,  as  before,  the  smallest 
angle  of  the  mesh.  If  0  <  f  (as  in  Fig.  5.1),  then  a  =  0  and  it  is  simple  to  prove  that  for 
this  problem 

(5.1)  effr  >  2.62  sin"  . 

On  the  other  hand,  if  0  >  f,  the  smallest  angle  is  a  =  f  —  0  and  in  this  case 

(5.2)  effr  <  6.86  sin^^  . 

Since  /  and  g  are  constant  and  1*2  coincides  with  u,  then  (4.9)  gives  for  this  problem: 

(5.3)  1.47  sin^^  <  eff  <  5.84  sin“^^j  . 

The  effectivity  indexes  (5.1)  and  (5.2)  corresponding  to  different  meshes  show  the  sharpness 
of  the  bounds  in  (5.3)  and  the  optimality  of  the  terms  sin^  (^)  for  their  dependence  on 
the  regularity  of  the  mesh. 
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6.  The  elasticity  problem.  We  shall  show  how  the  techniques  described  above  can 
be  applied  to  a  different  problem.  Let  us  consider  the  2D  linear  elastic  equations;  let  17,  Fn, 
Fa,  n,  Th,  Tj,  T  and  £  be  as  in  Sections  2  and  3;  let  Hfd(f7)  :=  {v  €  Hl(Q)2  :  v|Fd  =  0} 
be  the  space  of  admisssible  displacements;  let  e  and  a  :  H1(Q.)2  — >  R2*2  be  the  strain  and 
stress  tensors  defined  by: 


and 


.  .  1  /  dvi  dvj  \  .  . 

:”<v) 2  \dTj  +  dTJ  ' 


cr0(v)  :=  +2 fi£ij(v)  ,  i,j  =  1,2  , 


k=l 

where  A  and  n  are  the  Lame  coefficients  that  depend  on  the  Young  modulus  E  and  the 
Poisson’s  ratio  v  of  the  material: 

Eu  E  1 

A:=(T+„)(l-2„)’  ',:=20T7)’  •B>0’  °<l,<2- 


Given  a  body  force  f  €  L2(S 7)2  and  a  prescribed  traction  g  €  Zt2(rn)2  with  components 
locally  smooth  as  described  in  Theorem  3.2,  let  u  be  the  solution  of  the  boundary  value 
problem: 

{—(A  +  /i)V(divu)  —  fiAu  —  f  ,  in  17  , 
u  =  0  ,  on  Fa  , 

<r(u)n  =  g  ,  on  T„  ; 

For  v  and  w  €  Hpd(f7)  let 

a(v,w)  :=  /  V  <Tij(y)£ij( w)  ; 

a  is  a  continuous  symmetric  bilinear  form.  By  using  Korn’s  inequality  (for 
[20]),  it  is  proved  that  a  is  coercive  and  so,  the  energy  norm  ||  •  |jn  :=  a(-,  •) 7 
to  the  usual  Sobolev  norm  ||  •  ||j  n  on  Hpd(f7).  Problem  (6.1)  has  a  unique 
Hpd(17)  and  it  satisfies  the  weak  formulation  of  this  problem: 

(6.2)  a(u,v)=  /  f-v  +  /  g-v,  Vv€Hfd(S7). 

Jn  J  r„ 


instance,  see 
is  equivalent 
solution  u  € 


Let  ua  €  VA  :=  {v  €  Hpd(f7) :  \\T  €  V\{T)2,  VT  €  Th]  be  the  piecewise  linear  finite 
element  approximate  solution  of  problem  (6.2).  Proceeding  as  in  Section  2,  it  is  proved 
that  the  error  e  :=  u  —  u*  satisfies  the  residual  equation: 

«(e,v)=  /  f-v+  /  [g-<r(u*)n]  •  v+  V  f  [tr( u)nj,  ■  v  ,  Vv  €  (f7)  . 

J «  Jr-  ter-, Jt 
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For  any  triangle  T  €  Th  and  for  any  edge  t  €  Tn,  let  IITf  and  II,  g  be  the  local 
projections  of  the  data  defined  as  before;  let 

'Hu)n]e,  if/CTi, 

J/:=<  2{II<g-[a(ufc)n]|/}  ,  if^crn, 

,  0  ,  if  t  C  Ta  , 


and  let 


VT 


m2|n,f|2  +  i  J2  I£I2Ij'I! 


t€E- 


The  proofs  of  the  theorems  in  Section  3  can  be  immediately  extended  to  this  problem. 
Let  Up  €  {v  €  :  v|T  €  VP(T)2,  VT  €  Th)  be  the  approximate  finite  element  solu¬ 

tion  of  problem  (6.2)  in  this  space  and,  for  any  U  C  ft,  let  II  •  lit  :=  lu  £L»i  <’y(')£o(')- 

THEOREM  6.1.  With  the  definitions  and  assumptions  introduced  above 

£(CJ)%’|  +1IU-U 


Nils.  < 


(6.3) 


where 


Ten 
+  C 


pun 


Ter  /  \tec 


(6.4)  (CJ)2  := 


sup 

v€7>;  :  ||v||T*0 


]TT  /t(v  - y/)|  +lEteET  lm//(v~v/) 


□ 


THEOREM  6.2.  Let  w  €  Hpd(ft)  be  such  that  for  all  the  triangles  T  €  Th, 


and 


/(nrf)-w  =  |if|nrf|2 , 

Y  fJtw=Y  i£i2j2 

teET  Jt  teET 

3c;  >  0  :  ||w||r  <  CTnr  , 
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where  C'T  may  depend  on  the  shape  of  the  triangle  but  not  on  its  size  hT.  Then 


Vn  < 


(  sup  C 
\T€T* 


l  \ 

y 

Hello  +  C 

#  EKr 

+**  EMM 

\t€T  / 

\t€C  ) 

i 

□ 

The  constants  C£  and  C^,  can  be  computed  by  techniques  analogous  to  those  in 
Section  4.  They  depend  very  weakly  on  the  Poisson’s  ratio  v.  The  values  of  C'T  are  almost 
proportional  to  sin-^^*-)  (aT  the  minimum  angle  of  T )  as  for  the  Laplace  equation. 
Instead,  for  any  fixed  degree  p  >  2,  our  computations  show  that  C£  are  almost  proportional 
to  sin- 2  (“If-);  the  exponent  —  |  indicates  a  much  stronger  dependence  on  the  regularity 
of  the  mesh. 

Remark  6.1.  The  increase  of  the  factor  sin-1  (-f")  to  sin-*  (^-)  is  due  to  the 
constant  in  Korn’s  inequality.  Let  us  show  it  in  the  case  that  T  is  a  triangle  as  that  in 
Figure  6.1. 


C 


In  [22]  it  is  shown  that,  for  any  function  v  6  Hl(T)2, 

(6.5)  |v|;  r<Ci(^)2log(i^)  ||v|*  +  c,(^)Vi;,e , 

where  hT  is  the  diameter  of  T,  Q  is  the  biggest  circle  contained  in  T  and  r  is  the  length 
of  its  radius  (see  Fig  6.1).  The  estimate  (6.5)  is  optimal. 

If  were  used  in  the  denominator  of  (6.4)  instead  of  ||v||^,  the  term  sin-^  (^L) 

would  appear  as  in  (4.3).  On  the  other  hand,  for  functions  v  €  VP(T)2  with  three  degrees 
of  freedom  fixed  at  the  vertexes  B,C  (to  avoid  rigid  motions), 


II 


T 


< 


II 


V|| 


Q 


<  C3 


1 
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where  Q  is  the  quadrilateral  of  vertexes  B,  C ,  D,  E  in  Figure  6.1;  the  constant  C3  depends 
on  the  regularity  of  Q  (i.e.:  on  the  quotient  );  but  not  on  the  small  angle  aT.  Since 

for  aT  small, 


then,  for  these  functions  we  have  that  ^||'r  is  bounded  by  sin-1(^-)  (neglecting  in  (6.5) 
the  logarithmic  term).  Therefore,  since  in  (6.4)  the  supremum  can  be  taken  over  these 
functions  v,  we  can  expect  C'T  to  be  proportional  to  sin-?  (~2~)-  0 


The  following  table  gives  the  values  of  the  constants  C'Q  and  C£  in  the  estimate 


cavn-  C 


h*  (  E  lfl Z,T  )  +h*  E^U 

rer  /  \eec 


r 


<  Net 


n 


<  ca  *7n  +  !l«-«plln  +  c 


n 

2 


^2(Ei fi?,r  EigiL 

.Ter  /  \tec 


in  terms  of  the  minimum  angle  a,  for  p  =  2  and  for  different  values  of  the  Poisson’s  ratio. 


v  = 

.15 

v  = 

.30 

v  = 

.45 

a 

C'« 

cs 

C'Q 

C'Q 

31.34 

0.038 

30.31 

■ 

1 

11.10 

0.054 

10.81 

1 

0.075 

6.05 

0.066 

5.96 

1 

5.74 

0.088 

3.93 

0.078 

3.92 

1 

3.85 

37.5° 

0.100 

2.81 

0.089 

2.84 

0.054 

2.85 

45.0° 

0.113 

2.12 

0.101 

2.19 

0.061 

2.23 

52.5° 

0.125 

1.67 

0.113 

1.75 

0.069 

1.82 

60.0° 

0.136 

1.36 

0.124 

1.44 

0.077 

1.52 

Table  6.1.  Constants  of  equivalence  for  different  Poisson’s  ratios. 

7.  Conclusions  and  computational  aspects.. 

1.  The  error  estimator  can  either  underestimate  or  overestimate  the  true  error.  If  the 
solution  is  unsmooth  the  accuracy  of  the  estimator  could  deteriorate  (but  not  drastically — 
we  have  to  consider  a  higher  degree  p  in  (3.12)  and  the  deterioration  is  logarithmic) 

2.  The  main  factor  in  the  accuracy  of  the  estimator  is  the  geometry  of  the  elements. 

The  geometry  (angle  a)  has  to  be  understood  in  conection  with  the  differential  equation. 
For  example  when  an  elliptic  differential  operator  j=  1  2  a»>  dz  dx  (a0  constants)  is 
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considered,  the  equation  can  be  transformed  into  the  Laplace  equation  by  an  affine  trans¬ 
formation  which  will  modify  the  angles  of  the  triangles.  The  constants  arising  in  this  case 
are  those  of  the  transformed  mesh. 

3.  The  accuracy  of  the  estimator  depends  on  the  relation  of  the  axes  of  anisotropy 
of  the  solution  (i.e.:  the  eigenvectors  of  its  Hessian  matrix)  and  the  orientation  of  the 
triangles.  If  the  main  axe  and  the  orientation  of  the  triangles  are  orthogonal,  the  error  is 
overestimated;  instead,  it  is  underestimated  if  they  axe  parallel. 

4.  The  estimates  we  derived  are  theoretical  and  they  allow  us  to  define  correction 
factors;  for  example,  for  the  Laplace  equation  and  a  uniform  mesh  of  equilateral  triangles 
we  can  use  (from  table  4.1)  ^(0.121  •  0.850)  «s  0.32.  If  we  rather  needed  a  safe  estimator 
we  should  use  a  greater  corrector  factor  (say  1.5). 

5.  For  the  elasticity  equations,  our  estimates  show  a  larger  sensitivity  with  respect  to 
the  minimal  angle.  This  effect  grows  for  larger  Poisson’s  ratio. 

6.  In  practice,  the  bounds  on  the  effectivity  index  are  expected  to  be  better  than  in 
our  theoretical  analysis.  However,  (5.1)  and  (5.2)  show  that  they  cannot  be  much  better 
without  additional  restrictions.  Of  course,  the  examples  yielding  (5.1)  and  (5.2)  are  more 
or  less  extreme  cases.  For  a  detailed  computational  analysis  we  refer  to  [4], 
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